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Abstract. We investigate polyelectrolyte brushes in the osmotic regime using both theoretical analysis and 
molecular dynamics simulation techniques. In the simulations at moderate Bjerrum length, we observe that 
the brush height varies weakly with grafting density, in contrast to the accepted scaling law, which predicts 
a brush thickness independent of the grafting density. We show that such behavior can be explained by 
considering lateral electrostatic effects (within the non-linear Poisson-Boltzmann theory) combined with 
the coupling between lateral and longitudinal degrees of freedom due to the conserved polymer volume 
(which are neglected in scaling arguments). We also take the non- linear elasticity of polyelectrolyte chains 
into consideration, which makes significant effects as chains are almost fully stretched in the osmotic regime. 
It is shown that all these factors lead to a non-monotonic behavior for the brush height as a function of 
the grafting density. At large grafting densities, the brush height increases with increasing the grafting 
density due to the volume constraint. At small grafting densities, we obtain a re-stretching of the chains 
for decreasing grafting density, which is caused by lateral electrostatic contributions and the counterion- 
condensation process at polyelectrolyte chains. These results are obtained assuming all counterions to be 
trapped within the brush, which is valid for sufficiently long chains of large charge fraction. 

PACS. 61.25.Hq Macromolecular and polymer solutions; polymer melts; swelling - 36.20.-r Macro- 
molecules and polymer molecules - 61.20.Qg Structure of associated liquids: electrolytes, molten salts, 
etc. 



1 Introduction 

In recent years, polyelectrolyte brushes have been subject 
of extensive investigations both theoretically PEliniElEl 
IHITrTKniMlW and experimentally [TTHT31 

im2()U21U22U28U24U25U2n1l77] . Polyelectrolyte brushes are 
layers of charged polymer chains densely end-grafted onto 
surfaces of various geometries. They have important tech- 
nological applications notably in stabilization of colloidal 
dispersions in polar media, where forces between colloidal 
particles (coated by grafted polyelectrolytes) are modified 
and may be controlled in order to prevent fiocculation |281 
EU. Recent investigations piOll^lllllSllM7ll51l51ll0mmi2m3l 
ITlll5lll6lll7lll8llli)ll^0ll^lll^^3IIMII^5ll^6ll^7| have revealed 
the detailed phase behavior of charged brushes over a wide 
range of system parameters, namely, grafting density and 
charge fraction of chains. (Here we consider only planar 
brushes with no added salt.) It is known from scaling the- 
ories as well as from self-consistent field studies that the 
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interplay between steric, entropic and electrostatic contri- 
butions results in a variety of different scaling relations for 
the equilibrium height of the brush. 

In strongly charged brushes with a large effective charge 
density of grafted polyelectrolytes, i.e. when both graft- 
ing density and charge fraction of polyelectrolyte chains 
are sufficiently large, effectively all of the counterions are 
trapped inside the brush. Over a certain range of charge 
fractions and grafting densities, the repulsive osmotic pres- 
sure of counterions is the major effect that tends to swell 
the chains and balances their elastic stretching pressure. 
This regime is known as the osmotic brush regime, where 
the equilibrium thickness of the brush was obtained by 
scaling arguments to be independent of the grafting den- 
sity 3 4. The osmotic brush height may, however, exhibit 
a weak dependence on the grafting density as a result 
of non-uniform counterion distribution inside the brush. 
Specifically, variation of the counterion density profile in 
the direction normal to the anchoring plane (which may 
be caused by partial diffusion of counterions outside the 
brush layer) leads to a logarithmic dependence on the 
grafting density A weak dependence of the brush 

height on grafting density may also be produced by lat- 
eral inhomogeneities as will be addressed in this paper. 
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The crossover between these two mechanisms is regulated 
by the chain charge fraction and grafting density. In the 
regime of large grafting densities, steric effects dominate 
over electrostatic interactions and the excluded-volume 
repulsion (in good solvent conditions) balances the elas- 
tic pressure, leading to the so-called quasi-neutral brush 
regime @K]. The effects of excluded-volume interactions 
have also been studied in poor solvent conditions [5]. On 
the other hand, for larger charge fractions, the electro- 
static correlations due to Debye-Hiickel interactions be- 
tween mobile counterions and oppositely charged monomers 
may produce a dominant attractive pressure balancing the 
excluded-volume repulsion @K]. This leads to a collapsed 
brush regime, which has also been obtained in previous 
molecular dynamics simulations |T3tl4ll5j . 

In weakly charged brushes (which are not considered 
in this paper), the counterion cloud extends far beyond 
the brush height and the osmotic pressure of the counte- 
rions becomes irrelevant against uncompensated electro- 
static repulsion between charged monomers. This results 
in the charged or Pincus brush regime, where the elastic 
and the electrostatic pressures are balanced jB a ■ Details of 
the phase diagram and transitions between these regimes 
have been discussed extensively in the literature |T][2 ,3,4 

In recent molecular dynamics simulations |16| . struc- 
ture and equilibrium properties of salt-free planar osmotic 
brushes have been studied at moderate Bjerrum length. 
For moderate values of charge fraction and grafting den- 
sity, it has been observed that the counterions are mostly 
confined within the brush and generate an almost step- like 
density profile in the direction normal to the anchoring 
plane. In contrast, counterions exhibit a non-uniform dis- 
tribution in lateral directions. In this situation, the simu- 
lated brush thickness varies weakly with the grafting den- 
sity in contrast to the standard scaling law [HIE]- 

In scaling studies, it is commonly assumed that the 
counterions are distributed uniformly in lateral directions 
parallel to the anchoring plane and that the elasticity 
of chains is obtained from a linear (Gaussian) model. In 
this paper, we present a mean-field model for the osmotic 
brush, in which we consider corrections to both of these 
assumptions: Firstly, we take into account lateral electro- 
static effects including lateral variation of the counterion 
density profile around polyelectrolyte chains within the 
framework of the non-linear Poisson-Boltzmann (PB) the- 
ory, and secondly, we use a freely-jointed-chain model to 
mimic non-linear elastic stretching of the polyelectrolyte 
chains. Moreover, we take into account excluded-volume 
interactions between particles as well as the conserved 
polymer volume as the chain dimensions change, as will be 
described below. The results of our model for the behav- 
ior of the brush thickness, which we refer to as the non- 
linear osmotic brush regime, display a non-monotonic de- 
pendence on the grafting density. In particular, at moder- 
ate grafting densities, we find that the brush height weakly 
increases with the grafting density arising as a result of 
the interplay between lateral electrostatics and the cou- 
pling effects. This behavior will be compared with simula- 



tion results displaying a reasonable agreement. A weak de- 
pendence of the brush height on grafting density has also 
been observed in recent experiments |26II27| and has been 
compared with the non-linear osmotic brush predictions 
obtained using simple scaling arguments 120] • At small 
grafting densities, our model predicts re-stretching of the 
chains by lowering the grafting density, which is caused by 
an increasing electrostatic pressure acting on the chains in 
this limit. As we shall demonstrate, this behavior is regu- 
lated by the counterion-condensation process around the 
chains, which is captured within the non-linear PB the- 
ory. This regime has not been investigated in simulations 
and experiments yet, and we propose that one needs long 
chains to obtain this behavior (see the Discussion). 

An inhomogeneous distribution of counterions in lat- 
eral directions has been recently observed both in simu- 
lation [Tlll5lll6| and experiment |2Ull22U24ll25| . Analyzing 
experimental data on spherical charged brushes, Muller 
et al. |22j observe that the counterion distribution around 
polyelectrolyte chains indeed follows the the non-linear PB 
predictions as obtained within the cylindrical-ccll-model 
approach |.'i()U.'il j . Other experiments indicate that a large 
fraction of counterions binds strongly to the polyelectrolyte 
chains, thus counterions are not evenly distributed inside 
the brush 20,24,25 . On the other hand, the non-linear 
elasticity, which accounts for the finite extensibility of 
the chains, appears to be essential for strongly charged 
brushes at moderate grafting densities, as polyelectrolyte 
chains in these situations are found to be stretched up to 
60-80% of their contour length |16U18lliyil20ll22ll24ll26ll27| . 
Non-linear elasticity models have also been used in pre- 
vious works JU| (but without considering lateral effects 
as done here). It is important to note that the non- linear 
elasticity by itself can not generate a grafting-density de- 
pendence for the brush height |12| . 

We shall neglect the variation of the counterionic den- 
sity profile in the direction normal to the anchoring plane 
and assume that all counterions are trapped within the 
brush. In the range of parameters, where a sizable fraction 
of counterions could leave the brush layer, the step form 
of the density profile in the normal direction is changed 
and may be calculated using self-consistent field (SCF) 
techniques 1,2,6,8,10,11 . The SCF study of Zhulina and 
Borisov indeed shows that at large grafting densi- 
ties (osmotic regime), the variation of the density pro- 
file in the normal direction gives rise to a logarithmic 
dependence of the brush height on the grafting density. 
Such a dependence is obtained by treating the electro- 
static effects at the non-linear level using the PB equation 
but assuming that the concentration of counterions (and 
that of monomers) is smeared in lateral directions 
Our model, therefore, describes the complementary limit, 
which can be achieved by taking sufficiently long chains 
of large charge fraction (even at low grafting densities) as 
will be examined later. 

In the present study, we account for the volume inter- 
action between counterions and polyelectrolyte chains us- 
ing a closest-approach distance, but neglect the excluded- 
volume interactions between counterions themselves, which 
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may be justified in the considered range of parameters. In 
addition, the effective volume of the polyelectrolyte chain, 
which is not accessible for counterions, is assumed to be 
constant: when chains shrink (at a fixed grafting density), 
the available volume for counterions decreases leading to 
an increasing osmotic pressure. Such a coupling between 
longitudinal and lateral degrees of freedom is a very sim- 
ple way to mimic the back influence of the conformational 
changes of the chains on the true volume available for 
counterions inside the brush and their osmotic pressure. 

2 The geometry of the model 

The model which shall be employed in the present study 
to calculate the electrostatic free energy of the brush is 
based on the commonly-used cell model for rod-like poly- 
electrolyte solutions 32 33] - Each polyelectrolyte chain 
in the brush is assumed to have N spherical monomers 
of the same diameter bo among which a fraction of / is 
charged with a charge valency of q. We assume that each 
single chain is symmetrically enclosed in a cylindrical unit 
cell with radius D, which is determined by the grafting 
density p a . The chains length, L, is assumed to be much 
larger than the radius of the cell, L ^> D. Furthermore, 
each polyelectrolyte chain is modeled as a cylindrical rod 
(which we may refer to as the polyelectrolyte rod) with 
the same length and some radius R < D. We suppose that 
the total volume of the rod is constant, i.e., R 2 L = RqLq, 
where Lq and Ro stand for the length and radius asso- 
ciated with its fully stretched conformation; Lq is taken 
equal to the contour length of the chain, Lq = Nbg, but Rq 
can be generally different from the radius of the monomers 
due to different possibilities to choose a cylindrical rod 
model for a polyelectrolyte chain [HI] . We restrict our dis- 
cussion to a model in which Rq = bo/2 (Figure QJi). 

The electric charge of the polyelectrolyte chain, qfN, 
is assumed to be uniformly distributed over the surface of 
the rod with a linear charge density of r (all in units of el- 
ementary charge e). There is a fixed number of oppositely 
charged counterions, N c , with charge valency of q c and 
radius of r c confined inside each unit cell, such that the 
electroneutrality condition, tL = qfN = q c N c , is fulfilled. 
(The charge valencies q and q c , and also r, are defined to 
be positive.) Solvent is treated as a continuum background 
characterized completely by its dielectric constant e, and 
no additional salt is present in the system. 

In the simulation model, grafting points form a square 
lattice with the lattice spacing of 2D S i m and thus, a graft- 
ing density of p a = l/(4Z?^ m ). A natural space-filling 
choice of unit cell for this system is a square unit cell 
with the edge size of 2D S i m . Since the analytical solu- 
tion for the counterion distribution is only available for a 
cylindrical unit cell, we have to map the square simulation 
box to a cylindrical box. (Note that experimentally, the 
grafting-point distribution is highly irregular and thus dif- 
ferent from both simulation and analytical model.) There 
are different ways to adopt a cylindrical unit cell for such 
a planar brush as discussed in detail in Ref. Here we 
choose a cylindrical cell with the diameter, 2D, equal to 



(a) (b) 

\-2R -[ 




Fig. 1. The model used for the polyelectrolyte chain (a) and 
the cylindrical unit cell (b) for the brush as discussed in the 
text. The unit cell boundary in the simulations is shown by a 
square. 

the lattice spacing of the square lattice as shown in Figure 
[TJd, where D = D sim . 

3 Non-linear osmotic brush regime within 
mean-field approximation 

3.1 The electrostatic free energy 

The so-called mean-field or Poisson-Boltzmann theory of 
a cylindrical cell model for polyelectrolyte solutions has 
been presented long time ago 30,31,35 . In the mean- 
field approximation, the electrostatic correlations between 
neighboring cells are entirely neglected and the study of 
the system is therefore reduced to a single-cell study. In 
addition to this, correlations between counterions present 
in the same cell are systematically neglected, however, 
they remain still correlated with the polyelectrolyte rod. 
The system is then studied in its "ground state" , where all 
the fluctuations are neglected as well. Then, by virtue of 
the electroneutrality condition, it follows that the electric 
field vanishes over the cell boundary. 

A canonical field theory may be written for the system 
and the Poisson-Boltzmann (PB) theory is subsequently 
obtained as a saddle-point approximation to the corre- 
sponding action The canonical PB free energy per 
unit cell is calculated and can be written (in units of fcflT) 
as 

^f=- S ^ i /*[i(V^ + |5(r- W (,)] 
- A^ c ln[ / dv fl(r)e-*M] + C, (1) 
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where £b — e 2 /(AireeokBT) is the Bjerrum length (as- 
sociated with a medium with dielectric constant of e), 
£ = ^bQcT is the Manning parameter, and C is a constant, 
which contains contributions due to self-energies, kinetic 
energy of counterions and other numerical constants. The 
geometry function J7(r) takes into account the presence 
of hard walls and restricts the positions of mobile coun- 
terions to the cylindrical region of R < r < D, that is 
J?(r) = 1 for R < r < D and zero otherwise, where r de- 
notes the radial distance from the axis of the cylindrical 
unit cell. As a result of the saddle-point approximation, 
ST^f /Sip — 0, the dimensionless potential field i/;(r) is 
obtained to fulfill the so-called Poisson-Boltzmann (PB) 
equation 



V 2 i/> = ^5(r-R)-n 2 f2(r) 
R 



-V(r) 



(2) 



where k 2 = Anisqc^c/ f drf2(r) exp(— ip) is an unspeci- 
fied factor the value of which will be fixed once we specify 
the reference point of the potential 37 . To calculate the 
electrostatic free energy, one needs to solve Eq. (|2J for the 
potential field ip. (Here the physical electrostatic potential 
reads ksTtp / 'q c e.) Assuming that the polyelectrolyte rod 
is sufficiently long and that the solution for the potential 
field also has cylindrical symmetry, it follows that 



r^—\ — 2£ and f r— 

dr Jr=R. V dr /r=D 



0. 



(3) 



where we have used the global electroneutrality condition 

£L = N c £ B q 2 . (4) 

It is easily observed that the only non-trivial contribution 
to the PB free energy comes from the potential field in 
the interior region R < r < D, where the solution to Eqs. 
J21 and is available due to early works by Alfrey et 
al. [HOj and Fuoss et al. [21] ■ It was shown that the so- 
lution for ip(r) takes different functional forms depending 
on whether £ lies below or above a threshold value £ c . For 
a polyelectrolyte rod with fixed radius, R, the threshold is 
simply given by 



\n{D/R) 

T+Hd/R)' 



(5) 



while in our model with constant volume constraint, as we 
shall see later, £ c has to be determined from a transcen- 
dental equation. 

For £ < £ c , the solution reads 

2 2 

=ln[^sinh 2 (/31n^-tanh~ 1 /3)], (6) 

where (3 is given by the transcendental equation 

l-(3 2 



l-/?coth(-/31nf ) 



(7) 



(up to an additive constant and in units of k B T) as 

+ £] + ln[(£ - l) 2 - f3 2 } - \n(27T£ B fq c qNR 2 )}.(8) 
For £ > £ c , the solution to the PB equation is 

2 2 

4>(r) = H^r sin 2 (/31n ^ - tan" 1 /?)], (9) 

where (3 now satisfies the following equation 

1 + /3 2 



l-/?cot(-/31nf )' 



(10) 



Similarly, we can calculate the electrostatic free energy of 
the system, 

N q c X C ' ' R \ 1 + /3 2 

+ £] + Hit - l) 2 + P 2 ] - H^Bfq c qNR 2 )\ll) 

The free energies (|SJ and Ijlll) have also been derived by 
Lifson and Katchalsky _.35] using a charging process. 

We note that here R represents the actual radius of 
the polyelectrolyte rod and thus in our model is related to 
the Manning parameter of the system £ through 



(12) 



This dependence is induced by the volume constraint R 2 L = 
R 2 L , and the fact that the total charge of the rod is con- 
served, that means 



where 



£o = q c qf- 



(13) 



(14) 



is the Manning parameter associated with the fully stretched 
conformation of the rod. In the above formulation, coun- 
terions have been taken as point-like particles. To account 
for the finite size of counterions, we may define a clos- 
est approach distance between the rod and counterions, 
i? CQ , to be used instead of R in the preceding equations. 
Assuming that the counterions have radius of r c and re- 
calling the actual radius of the rod R from Eq. I|12|) . we 
have 

Note that £ is limited from above and below due to the 
geometrical constraints R ca < D and L < L respectively. 
As a result, one observes that 



£o < £ < Zu, 



Inserting the solution Eq. 10 into the expression we where 
obtain the electrostatic free energy per number of monomers 



2 /Rl 



(16) 



(17) 
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Finally, we remark that here the threshold Manning 
parameter, £ c , at which the functional form of the solution 
to the PB equation is changed, is determined from 



\nD/R ca {jc) 
l + lnD/R ca (^) : 



(18) 



in which R ca {^c) is given by Eq. ltT5|l . Equation l(T%|) is 
in fact obtained as an extension of Eq. JS} to the present 
model with constant volume constraint 1881 . 



3.2 The chain elastic free energy 

We assume a freely-jointed-chain (FJC) model to calculate 
elastic contributions to the total free energy of polyelec- 
trolyte chains in the brush. The exact free energy of such 
a chain model has a purely entropic origin and is obtained 
(in units of ksT) as 



T FJC , sinhy 

m h y coth y 



N 



1, 



(19) 



where N is the number of monomers and y is found from 

(20) 



L , 1 

— = coth y , 

La y 



in which L and Lq — Nbo are end-to-end distance and 
contour length of the chain respectively (see Appendix 
lA"|l . In terms of the modeled polyelectrolyte rod, L stands 
for the actual height of the rod and Lq for its height in 
the fully stretched situation. In the weak-stretching limit, 
Eqs. I|19|) and (|20|l lead to the free energy of a Gaussian 
chain 

TTFJC ot2 

(21) 



N 



2(Nb ) 



Whereas in the opposite limit of strong stretching, the 
model produces a non-linear elasticity with 



r "° -wi- L 



N 



Nbn 



const. 



(22) 



which is dominated by the finite length of the chain. It 
is clear that for almost-fully-stretched chains, deviations 
from the Gaussian behavior will be important |10| . 



3.3 Optimal brush height and its limiting behavior 

The total free energy of the brush per unit cell is the sum 
of the electrostatic and elastic free energies obtained in 
Eqs. © or lUTJ and 0, 



T 



T 



PB 



T 



FJC 



(23) 



The total free energy can be viewed as a function of the 
effective Manning parameter of the system, £, which varies 
according to the brush height L (see Eq. $l'6l ). It also de- 
pends on the system parameters, namely, charge fraction 
of the rod, /, valency of charged monomers, q, and that 




Fig. 2. Typical free energy of the present model (per number 
of monomers) is plotted as a function of the Manning parame- 
ter associated with the polyelectrolyte rod. Here it is calculated 
for / = 1, Ib = O.l&o, 9o = q = 1, r c = bo/2, and D = 1.56o- 
In this case, the lower and upper bounds for £ are £o = 0.1 
and £u = 0.4 respectively. The arrow shows the location of the 
optimal Manning parameter ~ 0.16. 



of counterions, q c , and finally on the ratio of the Bjerrum 
length, Ib, the cell radius, D, and the counterion size, 2r c , 
to the monomer size bo- 

The typical form of the total free energy (per unit cell) 
as a function of £ is shown in Figure where the free en- 
ergy is calculated for / = 1, Ib = O.l&o, q c = q = 1, 
i'c = bo/2, and D = l.bbo- As explained in Section l3~Tl the 
effective Manning parameter, £, is bounded from below by 
£o, and from above by given by Eq. I|17|). In this ex- 
ample, the lower and upper bounds on £ are respectively 
£o = 0.1 and £ u — 0.4. As seen in the graph, the total 
free energy has a minimum at an intermediate value of £. 
The reason is that by decreasing £ down to £cb the chain 
becomes highly stretched and its elastic free energy (see 
Eq. (|22Jl) and consequently its total free energy increases 
and eventually diverges due to its finite extensibility. In 
the other limit of large £, the available space for counte- 
rions decreases as £ — > (equivalently R ca —> D)\ thus 
the translational entropy of counterions dramatically de- 
creases leading to an increasing free energy in this limit. 
The optimal Manning parameter, £*, that minimizes the 
total free energy, lies somewhere between these two limits 
and corresponds to an optimal brush height L* = £o^o/£*- 

We thus obtain the optimal brush height, L*, by mini- 
mizing the total free energy of the system (I23|l for various 
grafting densities. The results will be compared with the 
simulation data in Section 0] In this Section, we focus on 
the generic predictions of the present model for the brush 
height and its limiting behavior. To this end and for the 
sake of simplicity, we consider a simpler version of the 
model, in which the counterions are taken as point-like 
particles (r c = 0). Hence, as in Eq. JTSJ) , 



R — R c 



(24) 
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Fig. 3. Log-log plot of the rescaled optimal height of the 
brush as a function of the grafting density for £b = O.lfeo, 
q c — q — 1 and point- like counterions (solid curves). Each of 
the curves corresponds to a charge fraction indicated on the 
graph. Dashed lines are the asymptotic behaviors obeying Eq. 
12811 for large grafting densities. 



We also take the charge valencies as q c = q = 1 . 

Figure |21 shows the optimal brush height plotted as 
a function of grafting density for the Bjerrum length of 
£b — O.lboj an d for three different values of charge frac- 
tion / = 1, 1/2 and 1/3 (solid curves). It is observed that 
the brush height has a non-monotonic behavior as a func- 
tion of grafting density: both in the limit of large grafting 
densities and small grafting densities, the brush height in- 
creases and eventually tends to its maximum value Lq — 
Nbo. Therefore, a lower bound is predicted for the equi- 
librium height of the brush that clearly depends on the 
charge fraction and the Bjerrum length. For charge frac- 
tions of / > 1/2 and is/bo ~ 0.1, this lower bound is 
about 50% of the contour length. This means that the 
polyelectrolyte rod remains increasingly stretched over a 
wide range of grafting densities. The limiting behavior of 
the brush height with grafting density can be understood 
both by asymptotic expansion of the free energy (see Ap- 
pendix[BJ) and by simple physical arguments as we present 
now. 

First, we consider the limit of large grafting densities 
Pa — > P™ ax > i n which p™ ax = 1/6q is the maximum graft- 
ing density. This corresponds to the limit of small cell ra- 
dius, i.e. D/Ro — > 1, where i?o = bo/2 in our discussion. 
In this limit, the available space for counterions inside the 
unit cell tends to zero and their entropy divergently de- 
creases. The produced osmotic pressure becomes the ma- 
jor repulsive pressure swelling the rod against stretching 
pressure of the elasticity, and bare electrostatic effects be- 
come negligible. The entropy of counterions is then well- 
approximated by the entropy of an ideal gas of particles, 
and can be written (per number of monomers and up to 
an additive constant) as 



where R is defined in Eq. 1)24(1 (see Appendix lB|l . In this 
limit, the chain has a large extension and its elastic free 
energy per number of monomers is given by Eq. I|22|l . 



felas 



ln(l-L/Z ). 



(26) 



Balancing the longitudinal pressures due to these two op- 
posing contributions, 



d_ 

dL 



[felas 



= o 



at fixed cell radius D, we find 



(27) 



(28) 



for pj)% - 
Section 



L 1 + / ' 

1, where, by the definition of the model in 



(29) 



The expression given by Eq. I(28|l is plotted in in Figure 
for the charge fractions / = 1, 1/2, 1/3 (dashed curves), 
where the coincidence with the predictions from minimiz- 
ing the full free energy in Eq. I(23|l (solid curves) is seen 
at large grafting densities. 

Note that within our model, the constant volume con- 
straint dominates when p a 6p — > 1, and the linear depen- 
dence on the grafting density in Eq. H28[) is induced by this 
constraint. It is important to note that the limiting behav- 
ior, Eq. 1(28(1 . identifies the behavior of the brush height 
if counterions around chains are taken as a uniform ideal 
gas of particles with no lateral electrostatic effects (see 
Appendix IB)) . But, as clearly seen from Figure El at lower 
grafting densities (e.g. at about p a &o ~ 0.1, which approx- 
imately corresponds to the simulated regime-see Figures 
and |H1 below), the predicted brush height (solid curves) 
deviates from the above limiting line (dashed curves) dis- 
playing a weaker dependence on the grafting density. This 
behavior at moderate grafting densities is induced by lat- 
eral electrostatic effects, which become increasingly im- 
portant and generate a minimum at intermediate grafting 
densities (Figure EJ. 

In the limit of small grafting densities p a ^o ^ 1) or 
equivalently D/Rq 3> 1, the present model is applicable 
for very long chains, since only in this case, counterions 
will be confined within the brush. In fact, to retain coun- 
terions inside the polyelectrolyte layer, one needs to con- 
sider chains of large charge fraction with the number of 
monomers N > (pa&o)" 1 ^ 2 ( see the Discussion). 

In this limit, the brush height shows different asymp- 
totic behavior in terms of the grafting density depending 
on whether the optimal Manning parameter, £*, is below 
or above the threshold value £ c as defined in Eq. |THJ). 
Thus, we shall distinguish two different scenarios. First 
we look at the case when the optimal Manning parame- 
ter of the rod is smaller than the threshold Manning pa- 
rameter < £ c . For our choice of parameters in Figure 



f]n[Tr(D 2 -R 2 )L], 



m 



= 0.1&o,<7 = q c = 1), this holds for p a bo < 0.1, 



(25) as can be checked from Eq. I|18l) . In such conditions, the 
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counterion cloud is highly diluted (counterions effectively 
de-condense |39I4()I4T| ) and there will be no electrostatic 
screening on the bare electrostatic potential of the rod 
ip(r) = 2£lnr [321 • This potential can be used to calculate 
the electrostatic energy (per number of monomers) 



Uele 



D 

R' 



(30) 



where R is defined in Eq. (|24|l and £ is related to the 
rod length L through Eq. H13|) . The residual entropic con- 
tribution of counterions may still be accounted for by 
assuming an ideal-gas entropy, s ci , as in Eq. The 
electrostatic free energy of the system is then written as 
Feiec/N Ri u e i ec - s cl which, for very large D/R , may be 
approximated by 



N 



/(£-2)ln 



D 
R' 



(31) 



This expression is indeed a leading-order term and can be 
derived by expanding the PB free energy © in terms of 
i?/D-powers (see Appendix |B|) . Using Eqs. l(3TJl and (|2"4~|) . 
the longitudinal electrostatic pressure can be calculated by 
differentiating T e iec with respect to L — £o^o/£- Since the 
polyelectrolyte rod is again highly stretched in this limit, 
it is reasonable to approximate the elastic contribution 
by the same non- linear (logarithmic) expression, Eq. l|2bl , 
as found in the strong-stretching limit. The equilibrium 
brush height is obtained by balancing these two contribu- 
tions to the longitudinal pressure on the polyelectrolyte 
rod (i.e. using Eqs. (|26|1 and J33I). The result (which has 
been calculated numerically) is shown in Figure^! (dashed 
curve) for / = 1 and £b = 0.1&o, where we also show the 
results from minimization of the full free energy, Eq. H23JI 
(solid curve a). The plot is made for p a down to 10 -6 6^~ 2 . 
For vanishingly small grafting densities p a 6§ — > 0, the en- 
tropic contributions become negligible compared with the 
bare electrostatic repulsion between monomers, and the 
equilibrium brush height behaves asymptotically as 



L /ln / 9 a 6g-2^ 



(32) 



This function is not shown in Figure 0] because it is valid 
for smaller grafting densities. 

In the second scenario, i.e. when the optimal Man- 
ning parameter associated with the chains becomes larger 
than the threshold > £ c , bare electrostatic interac- 
tions are partially screened as a result of the counterion- 
condensation process 39,40,41 . In fact, the PB electro- 
static potential (up to some logarithmic corrections) re- 
duces to the bare electrostatic potential of a rod with 
critical Manning parameter £m — T *-e., i/'( r ) = 21nr, 
when D /R — > oo |3*%ll4"2*] . The electrostatic energy of the 
system for D/Rq ^> 1 can be estimated using this poten- 
tial, which yields 



Uele 



V n R- 



(33) 



To estimate entropic contributions in this case, we may 
adopt the counterion-condensation picture .39 , 40 that only 




1e-06 



1e-04 



1e-02 



1e+00 



Pfio 



Fig. 4. Log-log plot of the rescaled optimal height of the 
brush as a function of grafting density. Solid curves show the 
results obtained from the minimization of the full free energy, 
Eq. 1231 . with / = 1, q c = q = 1, and point-like counterions, 
for a) l B = O.l&o (£o = 0.1), b) l B = 0.76 (fo = 0.7) and 
c) Ib = 1.2&o (Co — 1.2). The dashed and dot-dashed curves 
show the corresponding asymptotic estimates at small grafting 
densities for the cases a) (using Eq. and c) (Eq. 1551 ) 

respectively. The dotted line shows the brush height in the 
absence of lateral effects for / = 1 (Eq. 1361 '). 



a fraction of 1/^ of counterions are unbound and may con- 
tribute to the entropic pressure. Thus, the corresponding 
ideal-gas entropy of counterions, s C i in Eq. H25[). may be 
corrected by such a factor and used, together with Eq. 
B33JI . to derive the leading term of the electrostatic free 
energy, T e iec/N » u e iec - a«, as 



N 



l^R- 



(34) 



The expression l|34(l is confirmed again by a limit expan- 
sion of the PB free energy l|ll[) -see Appendix^] Calculat- 
ing the longitudinal electrostatic pressure from Eq. (|34|l 
and balancing it with the non-linear stretching pressure 
from Eq. H20I) . one finds that 



1 



2£o 



/In Pa fro' 



(35) 



The asymptotic expression, Eq. I|35|l . is shown in Figure 01 
(dot-dashed curve) along with the result from minimiza- 
tion of the full free energy, Eq. I)23|) (solid curve c) for 
a system with f — 1 and £b = l-2&o (Co = 1-2) for 
which the optimal Manning parameter, £„, remains al- 
ways above the threshold [35j. The above estimate, Eq. 
(|35fl . represents the result obtained from the model for 
Pa^o < At small grafting densities, the constant vol- 

ume constraint becomes unimportant. We also note that 
within our model, similar behavior is obtained for the 
brush height for the whole range of Bjerrum lengths, see 
the result for / = 1 and l B = 0.7b (f = 0.7) in FigureH 
(solid curve b). 
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The preceding discussion on the limiting behavior of 
the brush thickness demonstrates the important role of 
lateral electrostatic contributions to the total free energy 
that generate re-stretching of the chains at small grafting 
densities. Also we showed that the strength of these ef- 
fects is controlled by the counterion-condensation process, 
which is systematically included in the PB equation used 
in Sect ion |3~T1 33 41 42 . Both for weakly charged (£q < 1) 
and highly charged chains (£o > 1), lateral electrostatic 
contributions produce a repulsive longitudinal force acting 
on the chains, which logarithmically increases by decreas- 
ing the grafting density (Eqs. lf3"T)l and J23}). In particular, 
for highly charged chains, the force is independent of the 
brush height, i.e. —d!F e i ec /dL ~ In D/R, and results from 
the electrostatic screening due to condensation of counte- 
rfoils (see Eqs. (|34|l and {T3J). In any case, the increase of 
the brush height, which converges to the contour length, 
is logarithmically weak as the grafting density is lowered. 
We note that this behavior will not be obtained if lateral 
effects are neglected. In this case, the brush height remains 
independent of the grafting density, 



f 



1 + /' 



(36) 



which is shown by a dotted line in Figure 01 and follows 
from balancing the entropic term, Eq. (|25l) . and the non- 
linear elasticity, Eq. Ij26(l . where the volume constraint, 
Eq. (|23} , is also neglected 

As already mentioned, our model describes the situa- 
tion, in which counterions remain all inside the brush and 
distribute uniformly in the direction normal to the anchor- 
ing plane. In the range of parameters, when a considerable 
fraction of counterions leaves the brush (e.g. when chains 
are short or charge fraction is small), a different behavior 
will be obtained for the brush thickness as a result of the 
inhomogeneous distribution of counterions in the normal 
direction. In this case, if the concentration of counterions 
(and that of monomers) is assumed to be smeared out lat- 
erally, the brush height is obtained to decrease monoton- 
ically by lowering the grafting density as shown in Ref. 

Therefore, the behavior of the brush height in the 
presence of non-uniformities in both lateral and normal 
directions should be examined in an extended approach. 

Finally, we remark that the non-monotonic behavior 
of the brush thickness is not influenced by the elasticity 
model and qualitatively similar features are obtained when 
a Gaussian chain elasticity is used [4^j . 



4 Comparison to molecular dynamics 
simulations 

Computer simulations provide an excellent mean to study 
polymer systems. Extensive molecular dynamics simula- 




Fig. 5. A snapshot of the simulated polyelectrolyte brush with 
36 fully charged chains of N = 30 monomers (in light gray) 
at grafting density of p a = 0.12ct -2 . The Bjerrum length is 
l B = a. Counterions are shown by dark gray spheres. The box 
height perpendicular to the anchoring plane has been reduced 
for the sake of representation. 



In these simulations, a freely-jointed bead-chain model 
is adopted where the monomers are connected by non- 
linear springs with the so-called FENE (finite extensible 
non-linear elastic) potential, and end-grafted onto a rigid 
surface. The counterions are explicitly modeled as charged 
particles and no additional electrolyte is added. The sim- 
ulation box is periodic in lateral directions and finite in 
the z-direction normal to the anchoring surface at z = 0. 
We use the techniques introduced in Refs. QHJ^I to ac- 
count for the long-range nature of the Coulombic interac- 
tions in a laterally periodic system. The short-range re- 
pulsion between particles separated by distance r is mod- 
eled by a shifted Lennard-Jones (LJ) potential u(r) = 
4e L j{(<j/r) 12 - (cr/r) 6 + 1/4} for r/cr < 2 1 / 6 and u(r) = 
otherwise, with Lennard-Jones diameter of er being equal 
for both monomers and counterions. The counterions and 
charged monomers are univalent (q c = q = 1). For the 
simulations reported below, we choose moderate values for 
the Bjerrum length, i.e. £b = <J and £b ~ 2cr, which corre- 
spond to electrostatically-intermediate-coupling situation. 
With this choice of parameters, the average bond length 
(which is the result of the interplay between LJ repulsion 
and the FENE bond potential) is almost unaffected by the 
electrostatic repulsions and is bo = 0.98er. 

Figure [5] shows a snapshot from the simulation of a 



tions have been performed recently on polyelectrolyte brushesbrush with 36 chains of 30 monomers, which is fully charged 
at various grafting densities and charge fractions, both at and has a large grafting density of p a — 0.12a -2 . In this 
strong [13II14II15) and intermediate electrostatic cou- figure, the connectivity of the chains are preserved, so 
plings. that the chains may appear to extend beyond the sim- 
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Fig. 6. Density profiles of monomers p m (z) (open symbols) 
and counterions p C i(z) (filled symbols) as a function of the dis- 
tance from the anchoring surface. Shown are profiles for fully 
charged brushes of 36 chains of N = 30 monomers at is = <r 
and grafting densities (from bottom to top) p a o" 2 = 0.020 
(triangle-lefts), 0.042 (circles), 0.063 (squares), 0.094 (dia- 
monds), and 0.12 (triangle-ups) . 
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Fig. 7. Brush height as a function of grafting density for 
polyelectrolyte chains of N — 30 monomers (contour length 
Lq — 30a) with charge fraction of / = 1. Circles show the 
simulation data and squares are the predictions of the present 
mean-field model. The dotted lines (a) and (b) show the scal- 
ing predictions, Eqs. 1371 and 1361 . with Gaussian and non- 
linear elasticity respectively. Here Bjerrum length is is = o 
and charged particles are univalent. 



ulation box. (The simulation image has been produced 
using the VMD software package 29-) Simulated den- 
sity profile of monomers and counterions of the system 
in normal direction are shown in Figure El for the fully 
charged brush at several grafting densities. As seen, both 
monomers and counterions follow very similar nearly-step- 
like profiles with uniform amplitude inside the brush, which 
increases with grafting density (the monomers show a short- 
range ordering close to the anchoring plane, which is not 
relevant in the present study). These figures show that 
the counterions are mostly confined in the brush layer 
and that the electroneutrality condition is satisfied locally 
. One may observe (more clearly from simulated brush 
heights in Figures and |SJl that the polyelectrolyte chains 
are stretched up to about 60% of their contour length, 
and thus their elastic behavior is far beyond the linear 
regime. Therefore, within the chosen range of parameters, 
the simulated brush is in the strong-charging and strong- 
stretching limits and as we will show, it exhibits the non- 
linear osmotic brush regime. 

The average height of endpoints of the chains is one 
of the quantities which can be directly measured in the 
simulations. We compare the predictions of our theoreti- 
cal model for the brush height with the average height of 
the endpoints from our simulations in Figures {7\ and [H] for 
charge fractions / = 1 (with Is = a) and f — 1/3 (with 
£b ~ 2<t ) respectively. To obtain theoretical data from 
the model, we have taken into account the finite size of 
the counterions using r c — bo/2; moreover, we take bo = a 
while comparing these data with the simulation results. 

As seen from the Figures, the simulated brush height 
(circles) varies slowly with the grafting density, contrary to 
the predictions of scaling theories. According to the scaling 
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Fig. 8. Same as FigureQbut for / = 1/3 and Bjerrum length 
in « 2<7. 

studies 3,4,7], the osmotic brush thickness is obtained as 

Lose = -i=iV/ 1/2 6o, (37) 

which is independent of the grafting density (shown by 
dotted line (a) in the Figures). This relation is based on 
a balance between the osmotic pressure of an ideal gas of 
counterions inside the brush tt os = Nfp a /L, and the elas- 
ticity of Gaussian chains, ir e i as = 3p Sj L/Nb^ ) . If we adopt 
the non- linear elasticity for large stretching, Eq. I|26|l . the 
brush height is given by Eq. (|36|) . when lateral effects are 
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neglected. The result is also independent of the grafting 
density and shown by dotted line (b) in Figures and 
|H1 The present model considers corrections to the scaling 
theory by accounting for lateral electrostatic and coupling 
effects. The theoretical predictions (squares) appear to be 
in reasonable agreement with the simulation results. As 
shown in Ref. [34) by comparing different cell models for 
the brush, the deviations from the simulation data are in 
fact of the oder of the systematic error of the employed 
cylindrical cell model (Section EJ). In other words, choos- 
ing a different cell volume or a different rod model for 
the polymer chain leads to differences in theoretical pre- 
dictions that are of the same order as the deviations from 
the simulation data in Figures0and[S] Yet some of the dif- 
ferences may be in general due to excluded- volume effects 
that have been treated in a simplified way in our study. 
In particular, we have neglected the excluded-volume re- 
pulsion between counterions; a rough estimate of these ef- 
fects, however, shows that such volume effects are almost 
one order of magnitude smaller than the osmotic pressure 
of counterions (see the Discussion). 

Another effect which may play a role in the simulated 
systems is lateral wiggling of the polyelectrolyte chains. 
In fact, the volume constraint, that has been used in our 
model, can account for some of these effects in a very 
simple way, since the effective chain radius increases for 
decreasing brush height. However, conformational changes 
of the chains may affect the validity of the cell model that 
was employed for calculating the electrostatic contribu- 
tions to the free energy. From lateral monomer density 
profile obtained within the simulations j47) . we conclude 
that wandering of chains is well-confined to a cylindrical 
region, the size of which is smaller than the "cell radius" . 
The cell radius is close to the decay range of the lateral 
monomer density profile only for large grafting densities, 
thus in this regime, lateral wiggling of the chains might 
be responsible for some of the deviations between the the- 
oretical predictions and the simulation data. 



5 Conclusion and Discussion 

We have studied strongly charged polyelectrolyte brushes 
in the non-linear osmotic regime. In order to bring out 
the features of the non-linear osmotic brush regime more 
clearly, we chose a moderate Bjerrum length of the order 
of the monomer size in our simulations. Nonetheless, the 
present theory is obtained for the whole range of Bjerrum 
lengths, and similar features are obtained at large, more 
realistic values of the Bjerrum length (see Section EU- 

In the simulations PB]> we have observed that the 
brush height varies slowly with the grafting density of the 
brush. In this situation, polyelectrolyte chains are strongly 
stretched and counterions are mostly confined within the 
brush layer. In the direction normal to the anchoring plane, 
counterions produce an almost step-like density profile, 
while in lateral directions, an inhomogeneous distribu- 
tion is observed around polyelectrolyte chains at moderate 
grafting densities [14115116) . 



We have investigated this system using a theoretical 
model, in which we consider a number of corrections to 
the standard scaling theories. Firstly, the electrostatic free 
energy of the system is calculated within the non-linear 
Poisson-Boltzmann theory, which includes lateral electro- 
static effects. On the other hand, the elastic free energy of 
the chains is calculated using a freely-jointed-chain model 
providing a non-linear elasticity in the strong-stretching 
limit, which accounts for the finite extensibility of the 
brush chains. We have taken into account also the back in- 
fluence of the conformational changes of the chains on the 
osmotic pressure of counterions. Such a coupling between 
longitudinal and lateral degrees of freedom is modeled in a 
simple way by assuming a constant effective volume for the 
chains. The results of our model display a non-monotonic 
behavior for the brush height as a function of the graft- 
ing density. In particular, at moderate grafting densities, 
we find that the brush height weakly increases with graft- 
ing density, which arises as a result of the interplay be- 
tween lateral electrostatics and the coupling (volume con- 
straint) effect. This behavior agrees with the simulation 
data within the systematic error of our model; the devi- 
ations may originate mainly from non-uniform boundary 
conditions in the simulations. Predictions of our model 
also agree with recent experiments on osmotic brushes on 
a semi-quantitative level as discussed in Refs. [261127) . 

As already discussed, the behavior of the brush height 
at very large grafting densities (closely-packed brush) is 
dominated by the volume constraint used in our model. 
The characteristic of this regime is the linear increase of 
the brush height with grafting density (see Eq. 1)28)1 ~). In 
this regime, excluded-volume effects become relevant and 
should be considered explicitly in the balance equations. 
While, in our study, volume interactions are only effec- 
tively modeled. Nonetheless, at lower grafting densities 
(e.g. at about p a b^ ~ 0.1 as investigated in the simu- 
lations), lateral electrostatic effects involved in the PB 
free energy become increasingly important and lead to a 
weaker increase of the brush height with grafting density 
(see Figure OJ- 

For large grafting densities, possible contribution of 
the excluded-volume repulsions between counterions may 
be estimated as follows. Effectively, the second- virial con- 
tribution from counterions to the total free energy (per 
unit cell and in units of fcgT) may be written as 

?v2 1 _2 T/ 1 Nf 

— = ^V2 P^oell = ^\ {D ,_ R , )V (38) 

where V2 > is the effective virial coefficient (for good sol- 
vent condition) , and ~p ci is the average density of counteri- 
ons in a unit cell of volume V ce u — ir(D 2 — R 2 )L. The lon- 
gitudinal pressure coming from the excluded-volume re- 
pulsions between counterions, ity™ 9 , is calculated by dif- 
ferentiating Eq. 1)38)) with respect to L, the brush height. 
The produced pressure is then compared with the longi- 
tudinal osmotic pressure of counterions obtained from Eq. 
PS), that is 7T^ on9 ~ f/L. We find 
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where we have also considered the closest approach dis- 
tance between counterions and monomers as discussed in 
the text (here we have taken univalent particles for sim- 
plification). Now, for largest grafting densities used in the 
simulations, e.g. pj% ~ 0.1, and for fully charged chains 
/ = 1 (where L/L ~ 2/3), it follows that ir^ n9 /tt%™9 ~ 
(l/lO)^/^ C ; where v§ = 7r6o/6 is the hard-core sec- 
ond virial coefficient. This ratio is even smaller for smaller 
charge fractions and grafting densities, and roughly yields 
the largest estimate based on the simulation data. Clearly, 
if we assume that t>2 ~ t>2) the excluded- volume pres- 
sure of counterions, which is neglected in our model, is 
found to be almost one order of magnitude smaller than 
their osmotic pressure. 

At small grafting densities, our model predicts a weak 
re-stretching of the chains by lowering the grafting den- 
sity, which is produced by lateral electrostatic contribu- 
tions and is controlled by lateral rearrangement of counte- 
rions (countcrion-condensation process) around the chains 
fSection I3.3fl . These results are valid only if the optimal 
height of the brush is much larger than the distance be- 
tween neighboring grafted chains. (In particular, the coun- 
terion condensation occurs in the presence of long chains, 
when end effects are negligible.) This situation may be 
achieved by taking polyelectrolyte chains with sufficiently 
large contour length and at large charge fractions. Oth- 
erwise, the cylindrical symmetry, assumed in the calcu- 
lation of the electrostatic free energy, and also the os- 
motic conditions are no longer retained. We may estimate 
the regime of validity of our results as follows. At small 
grafting densities, a large fraction of counterions tends to 
leave the brush and consequently, the system is pushed 
away from the osmotic regime: for large charge fractions, 
it will enter the Pincus brush regime, where uncompen- 
sated electrostatic repulsion between monomers balances 
the elastic stretching, and for small charge fractions of the 
chains, mushroom conformations, that are entropically fa- 
vorable, are formed by shrinkage of the chains onto the 
anchoring plane. The crossover from osmotic to Pincus 
brush regime can be identified by introducing the average 
height of the counterion layer over the anchoring plane. 
The thickness of the counterion layer is in fact character- 
ized by the effective Gouy-Chapman length of the plane, 
that is Xqc = l/27r£_B./V/p a (we restrict our discussion to 
the case with univalent particles). In the osmotic brush 
regime, where the surface charge density Nfp a is large, 
one has Xgc < L, whereas in Pincus brush regime with 
small surface charge density Xqc 3> L. (In the simula- 
tions, we have Xqc/L ~ 10~ 3 for / ~ 1 and p a b\ ~ 0.1.) 
To retain counterions inside the brush, one needs to take 
polyelectrolyte chains with number of monomers N larger 
than a certain threshold number AT* . To estimate AT* , one 
can use scaling arguments that yield Pincus brush height 
as 

L PB ~ NHsfpJbl. (40) 

Now, given the scaling relations for the osmotic and Pincus 
brush regimes, Eqs. (|3"?| and l|4T))) respectively, one can 
find the boundary relation between these two regimes as 
j\j-2 f 3 / 2 l B p a bo ~ 1 The threshold number of monomers 



A* follows from this relation, i.e., 

n, = r"/v&S)- 1/a (£) 1/a - (4i) 

One notes that for / ~ 1 and £b ~ bo this corresponds 
to the overlapping threshold for neighboring chains, that 
is when p a L^ ~ 1 with Lq = Nbo as the contour length 
of the chains. To give a numerical estimate of N*, we use 
the simulation parameter £b — bo for fully charged chains 
(/ = 1): at grafting densities as small as p a b\ ~ 10~ 6 , Eq. 
(|4T)) gives A 7 , ~ 10 3 , while for pJ>o ~ 10~ 2 we would have 

w 10. Note that the latter case coincides with the min- 
imum grafting density used in the simulations (see Section 
EJ, where polyelectrolyte chains bear N — 30 monomers. 

Now, to prevent formation of mushrooms at small graft- 
ing densities and small charge fractions, one may take 
longer chains with N > AT,, = (yO a 6g) _1 , which is a more 
stringent condition on N, i.e. > A*. This corresponds 
to the condition that the Gaussian size of the polyelec- 
trolyte chains, A fl / 2 6o, be larger than the distance between 

— 1/2 

neighboring chains p a . (Note that we have assumed a 
Gaussian polymer, thus our estimate gives an upper bound 
for iV*«.) At small grafting density of psb\ ~ 10 -6 , one has 
Af,* ~ 10 6 , a quite large value, and for p a b\ ~ 10~ 2 , we 
obtain Af*» ~ 100. We conclude that the regime, where 
the brush height has been predicted to increase with de- 
creasing the grafting density is observable in experiments 
and simulations by choosing long enough chains. 

In this paper the system is studied within a mean- 
field approximation, which is valid only if electrostatic 
correlations are negligible. Such effects become important 
at strong electrostatic couplings 48 , i.e. at large Bjer- 
rum lengths and charge valencies, and can produce sig- 
nificant attractive pressure on the chains. In such condi- 
tions, a scaling theory has been developed which predicts 
a collapsed brush regime observed also in simulation ^21 • 
At moderate Bjerrum length and for univalent particles 
as chosen in our simulations, correlation effects are not 
yet important, and as we showed, the predictions of the 
present mean-field theory are very close to the simulation 
results. 

An interesting problem is to extend the present results 
to include the variation of the counterionic density profile 
in the direction normal to the anchoring plane, which may 
be caused by the diffusion of counterions outside the brush 
at low grafting densities or at small charge fractions. The 
self-consistent field analysis of Zhulina and Borisov 
reveals that if these effects are accounted for using the 
non-linear PB equation, the brush height monotonically 
decreases by lowering the grafting density provided that 
the concentration of counterions (and that of monomers) 
is assumed to be smeared out laterally. Our results show 
a re-increase of the brush height by lowering the graft- 
ing density, when the counterion profile is assumed to be 
uniform in normal direction, and allowed to admit a later- 
ally inhomogeneous form according to the non-linear PB 
equation. An extended approach should, therefore, exam- 
ine the interplay between these two mechanisms. 
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Finally, we remark that similar calculations can be 
done for slightly different models |34| . e.g. one may con- 
sider a volume charge distribution for polyelectrolyte chains 
Clearly, the Coulomb self-energy of the polymer would be 
larger in this case and a larger brush height is predicted 
as compared with the studied model. 



the free energy T FJC /N w 3L 2 /2(Nb ) 2 ■ In the strong- 
stretching limit boF 3> 1, a non- linear force-extension re- 
lation is reached from Eq. (fHjl . that is L/b N 1— 1/boF. 
Using this, we find the elastic free energy T FJC /N sa 
- ln(l - L/Nbo) + const. 
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B Asymptotic expansions of the PB free 
energy 

Here, we briefly discuss the asymptotic expansions of the 
PB free energy both in the limit of small and large grafting 
densities. For simplicity, we assume q c = q. 

B.l Small grafting densities (L>/i?> 1): 



A The free energy of a freely-jointed chain 

We calculate the free energy of a freely-jointed chain con- 
sisting of N links (monomers), each of fixed length bo and 
able to point in any direction independently of each other. 
It is convenient to perform the calculation in an isobaric 
ensemble, in which the chain is considered to be stretched 
by applying a constant force F (in units of kgT). Since 
the monomers are assumed as rigid links, the configu- 
ration space of the chain is spanned by a set of angles 
{9i,ipi} specifying orientations of the monomers labeled 
by i = 1,...,N (a standard spherical frame of coordi- 
nates is chosen with z-axis pointing in the same direction 
as the force). The partition function of such a polymer 
chain can be written as 



Z F = 



271 



o 2?r Jo 2 



dd 



■ sin v e 



Integrating Eq. I|42|) we find 



Z F = 



J> F 



-b F 



2b F 



N 



(42) 



(43) 



Now, the extension (end-to-end distance) of the chain, L, 
is calculated from L = d In Z F / dF as 



L 1 

= coth bnF — . 

b N b F 



(44) 



In this case, there are different limit expansions for the PB 
free energy depending on whether the Manning parameter 
is smaller or larger than £ c . Note that in this limit, the 
threshold Manning parameter itself tends to one, i.e. £ c — * 
1 (see Eq. |[TSJ}, also Ref. [3T])- 



B.l.l Case with £ < £ c . 

Inspecting Eq. J7J, one can find that in the limit of D/R — > 
co, the solution for (3 tends asymptotically to (3 = (£ — 1). 
Hence for D/R ^S> 1, we can propose the following form 
for (3, 

/3 2 «(£-l) 2 (l-x), (46) 

where a; is a small function of £ and D/R, and may be 
determined as follows. Rearranging Eq. Q, we can find 
a more convenient equation for the forthcoming limit ex- 
pansions, 

1 R 2 1 + /3 2 (£-i)-0 y ' 

Now using Eq. 146|) into Eq. (|47|l and expanding for small 
x, we have 



4£ ( D 



2-i\R 



2(«-l) 



(48) 



Similar limit expansions can be performed for the poten- 
tial field, Eq. ©, and the PB free energy, Eq. (JjJJ, with x 
given by Eq. JtHJ, which lead to 



Using Legendre transformation we can calculate the iso- 
choric free energy of the system (in units of fcgT), T F,1C — 
— \o.Zp + LF, hence 



T 



FJC 



N 



ln^M+fcoFcothfroF-l, (45) 
OqF 



which was used together with Eq. 144fl in Section 

In the weak-stretching or Gaussian-chain limit b$F <^ 
1, the end-to-end distance is obtained by a proper expan- 
sion of Eq. lB4^l asL~ Nb^F/S, which is then used to get 



^(r)«2£ln-, 
(where we have assumed ip(R) = 0), and 



T 



PB 



N 



/(£-2)ln^, 



(49) 



(50) 



as estimated in Eq. (|31(l in the text. Clearly, both bare 
electrostatic repulsions of the charged rod and entropy of 
mobile counterions contribute in Eq. (|5U|) respectively as 
f^hxD/R and —2f\nD/R. However, in the limit of very 
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large D / R, the longitudinal entropic pressure of counteri- 
ons becomes vanishingly small relative to the longitudinal 
bare electrostatic pressure. This can be seen from Eq. (|50|) . 
noting that the longitudinal pressure is obtained by dif- 
ferentiation of the corresponding term in the free energy 
with respect to the cell volume ttD 2 L at fixed cell radius 
D. Formally, the longitudinal osmotic pressure, 7T^ s °™ 9 , will 
be Tta° n9 = —dT e i ec /irD 2 dL, where L is related to £ via 
Eq. (JT3J) and T eUc is given by Eq. I|50l) . 

In contrast, both bare electrostatic repulsions and en- 
tropic effects produce lateral pressures of the same leading 
order. (The lateral pressure is obtained by differentiation 
of the corresponding term in the free energy with respect 
to irD 2 L at fixed rod length L.) The latter result is known 
as Manning limiting law for osmotic coefficient of dilute 
solutions of weakly charged polyelectrolytes [101 • This lim- 
iting law follows directly from Eq. I|50[) as we explain now. 

For rod- like polyelectrolytes with £ < £ c , the Manning 
limiting law states that the osmotic coefficient, ^, tends 
to a finite value as the solution becomes highly diluted, 



lim 

_D/_R-»c 



"=(l-f). 

) z 



(51) 



The osmotic coefficient here is defined as the ratio of the 
lateral osmotic pressure acting on the cell boundary, 



_Lat 



(D) 



2ttDL8D 



(52) 



to the lateral osmotic pressure of an ideal gas of particles, 
Pid, filling the cell under similar conditions, i.e. 



r Lat 



Pi 



id 



(53) 



Now, given the ideal-gas pressure Pid = N c /ttL(D 2 — R 2 ) 
and the (PB) electrostatic free energy of Eq. I|5()(l . one can 
easily recover the limiting law, Eq. (|51|l . using Eqs. H52[) 
and 



B.1.2 Case with £ > £ c . 

In this case, it follows from Eq. HIQfl that /3 tends to zero 
(as (In D/R) -1 ) when D/R — > oo. Using this, one can 
observe that the potential field of Eq. © reduces to 



^(r)w21n- + 21n[l + (£-l)ln-] 



R 

and the PB free energy, Eq. Qlljl- to 



R 1 



PB 



-.-(in* 
N £ R 



(54) 



(55) 



as estimated heuristically in Eq. i|34|) in the text. Clearly 
in this case when £ > £ c , unlike the case with £ < £ c , elec- 
trostatic repulsions and entropic effects have contributions 
of the same order of magnitude in the total longitudinal 
pressure which acts on the rod. This result may be under- 
stood in terms of the counterion condensation model 0U] 
as explained in Section 13.31 



The Manning limiting law for the (lateral) osmotic co- 
efficient is now obtained (from Eqs. I|52|l . I|53|l and (|55[l ) 

as 



lim v 



1 



(56) 



where £ > £ c . 

In our model where R is not fixed but depends on £ 
(see Eq. (|15[l ). the preceding discussions hold when D ^> 
r c + R Q . 



B.2 Large grafting densities (D/Rm 1): 

In a cell model with fixed R, the threshold Manning pa- 
rameter £ c tends to zero as D/R — > I [31]. Therefore, for 
finite values of the Manning parameter £, we have to use 
Eqs. {TDJ and ([TTJ for £ > £ c . 

Starting from Eq. 1)1 0|l. we can find an approximate 
expression for [3 in the limit of D/R — > 1. Defining 8 = 
D/R — 1 and expanding Eq. (|f Oj l for small 8, we obtain 



f 



(57) 



Now replacing (3 2 in Eq. (|1I|I and expanding in terms of 
8, we have (up to an additive constant independent of 8) 



T 



PB 



N 



-flnS + 0(S). 



(58) 



This is the entropic free energy of an ideal gas of particles 
(Eq. lHH ) up to the leading order, and shows that in this 
limit, the main contribution to the PB free energy comes 
from the entropy of counterions. 

In our model with the constant volume constraint, this 
limit (namely D — > r c + Rq) has to be handled with care. 
In fact, the upper limit on £, that is £„, tends to the lower 
limit £o 5 and so does the optimal Manning parameter. The 
threshold £ c becomes smaller than £o 5 therefore, the sys- 
tem indeed satisfies £ > £ c condition [SHj, so that the 
above discussion remains valid. 
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